Dominant superconducting fluctuations in the one-dimensional extended 

Holstein-extended Hubbard model 
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The search for realistic one-dimensional (ID) models that exhibit dominant superconducting (SC) 
fluctuations effects has a long history.—^ In these ID systems, the effects of commensurate band 
fillings-strongest at half-filling-and electronic repulsions typically lead to a finite charge gap and the 
favoring of insulating density wave ordering over superconductivity. Accordingly, recent proposals 5,6 
suggesting a gapless metallic state in the Holstein-Hubbard (HH) model, possibly superconducting, 
have generated considerable interest and controversy, with the most recent work demonstrating that 
the putative dominant superconducting state likely does not exist.— ~— In this paper we study a model 
with non-local electron-phonon interactions, in addition to electron-electron interactions, this model 
unambiguously possesses dominant superconducting fluctuations at half filling in a large region of 
parameter space. Using both the numerical multi-scale functional renormalization group for the full 
model and an analytic conventional renormalization group for a bosonized version of the model, 
we demonstrate the existence of dominant superconducting (SC) fluctuations. These dominant SC 
fluctuations arise because the spin-charge coupling at high energy is weakened by the non-local 
electron-phonon interaction and the charge gap is destroyed by the resultant suppression of the 
Umklapp process. The existence of the dominant SC pairing instability in this half-filled ID system 
suggests that non-local boson-mediated interactions may be important in the superconductivity 
observed in the organic superconductors. 

PACS numbers: 71.10.Fd, 71.30.+h, 71.45.Lr 



I. INTRODUCTION AND BACKGROUND 



Many interesting novel electronic materials, includ- 
ing charge-transfer solids and conducting polymers, are 
quasi-one-dimensional physical systems and exhibit clear 
effects of both electron-phonon (e-ph) and electron- 
electron (e-e) interactions. Accordingly, a large number 
of one dimensional (ID) Hamiltonians have been studied 
as microscopic models for these systems, with consid- 
erable success in describing the range of ground states, 
nonlinear excitations, and optical and transport proper- 
ties of these materials. One outstanding challenge for 
theorists has remained, however: namely, the identifica- 
tion of a ID model, incorporating both electron-phonon 
and electron-electron interactions with realistic interac- 
tion parameters that exhibits dominant superconduct- 
ing fluctuations. This has proven especially difficult be- 
cause of the strong tendency in these ID systems to form 
insulating density waves ground states: charge density 
wave (CDW), spin density wave (SDW), or bond-order 
wave (BOW). Hence the recent suggestion^ that the ID 
Holstein-Hubbard (HH) model at half-filling might have, 
for a narrow window of parameters, a region of dominant 
superconducting fluctuations has generated considerable 
interest. Although some later analyses^ - — suggest that 
this model likely does not have dominant superconduct- 
ing fluctuations, the initial results^ provided support for 
the existence of a novel gapless metallic phase intermedi- 
ate between the familiar "Peierls" (CDW insulating) and 
the "Mott" (SDW insulating) phases. This result, tanta- 
lizingly close to the long-sought dominant SC state, has 



rekindled interest in finding a ID model which does have 
clearly dominant SC fluctuations for a range of physically 
reasonable parameters. 

In this paper, we define and study such a model, 
which we call the "extended Holstein-extended Hubbard" 
(EHEH) model. We will define the model precisely in 
the ensuing section, but for now we note simply that it 
involves both on-site (U) and nearest-neighbor (V") re- 
pulsive e-e interactions and non-local e-ph interactions. 

Before describing the EHEH model in detail, let us first 
briefly summarize the results in the Holstein-extended 
Hubbard (HEH) model, as these serve as a very useful 
guide to the surprising subtleties encountered in studies 
of this seemingly simple model. The Hamiltonian for the 
HEH model is 

Hheh = -*5Z(cJ+i,<7 < V + H - c ) + w o X] a i ai 

i.a i 

i i 

+ ge P y^(a| +ai)ni, a . (1) 

i,a 

Here electrons move in a tight-binding ID lattice, where 
c l a ( c i.<y) creates (annihilates) an electron at site i with 
spin <7, riiff is the electron number operator, ni = n^-f + 
Tii,i, and aj creates a (dispersionless Einstein) phonon of 
frequency luq at site i. The electron density is coupled lo- 
cally with strength g ep to phonons on the same site. The 
e-e interactions are described by the standard extended 
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Hubbard model with on-site interaction U and nearest 
neighbor interaction V. 

To recall the results for this model we begin with the 
case g cp = 0, so the model reduces to the familiar ex- 
tended Hubbard model^ which incorporates only e-e in- 
teractions. We focus on the half-filled case in the phys- 
ical region with U and V both positive. For 2V >> U, 
CDW fluctuations are dominant and the ground state 
is indeed a long-range ordered CDW. In contrast, for 
2V << U , spin density wave fluctuations are dominant 
and the ground state is an (algebraically decaying) SDW. 
Recent results^— have shown that there also exists a 
small region around the line U = 2V in the U — V plane 
in which the ground state is a long-range ordered BOW 
(also sometimes called a bond-charge density wave). But 
there is no evidence for dominant superconducting fluc- 
tuations in this model for half-filling with U and V pos- 
itive. This result coincides with the naive expectation 
that for purely repulsive interactions, superconducting 
pairing should not exist. For the pure Hubbard model 
(V = 0)the exact solution^ shows that the charge gap 
opens immediately for U > but that it is exponentially 
small for small U. 

We consider next the case of e-ph interactions only, 
in which we recover the equally familiar ID Holstein 
model. 19 In a pioneering study nearly three decades ago, 
Hirsch and Fradkin 20 argued that for spin 1/2 electrons, 
in the half-filled Holstein model a charge gap opens un- 
conditionally for any non-zero e-ph coupling and at any 
finite phonon frequency. In the adiabatic limit (ujq — 
0) the gap is also exponentially small in g cp . While 
some subsequent studies have supported this results^ - — 
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have suggested that for sufficiently large luq 



a finite value of g cp is required to open the gap, so for 
sufficiently small g cp there is a gapless phase. 

Similarly, in the full HEH model, although the ini- 
tial study^S suggested no gapless phase, more recent 
studies^ ' 27 ' 28 ' 31 " — , have found the "gapless" metallic 
phase mentioned above in a small region of the g ep -U 
plane around the line U = 2g% p /ujo- Ref[6| gives a plau- 
sible intuitive interpretation of this phase, as well as a 
clear statement of the important caveat of the likely ex- 
ponentially small nature of any gap in this region makes 
drawing definitive conclusions from these numerical stud- 
ies very difficult. A further reflection of the subtleties in 
this problem is the fact that other studies, including some 
recent results, do not find this gapless phase^2r— 

Whatever the ultimate resolution of this issue in the 
HH model,— the "fragility" of the gapless metallic state 
and the absence of dominant superconducting fluctu- 
ations in that model strongly suggest that we seek a 
more robust model in which it is clear that there are 
dominant superconducting fluctuations. A hint as to 
what sort of model could produce this result comes from 
considering a problem simpler than the full many-body 
problem: namely, the two-electron "bipolaron" prob- 
lems. Recently, several groups have studied this rela- 
tively more manageable problem in models with long- 



range e-ph interactions^ 3 - - — One of the main differences 
between the Holstein model, which has only on-site e-ph 
interactions, and models with non-local or long-range e- 
ph interaction is that in the latter models the mass of the 
bipolaron is reduced considerably^ 3 - - — with the possible 
enhancement of superconductivity. If these results trans- 
late to the many-electron problem, then one should ex- 
pect non-local e-ph interactions to enhance the possibility 
of superconductivity and lead to dominant superconduct- 
ing fluctuations in a large region of parameter space. Ac- 
cordingly, in the ensuing sections, we study a model-the 
"extended Holstein-extended Hubbard" (EHEH) model- 
which includes both non-local e-ph interactions and the 
e-e interactions of the extended Hubbard model. We shall 
see that this model does indeed allow for dominant SC 
fluctuations in a substantial portion of the U — V plane. 



II. THE EHEH MODEL AND THE MFRG 
METHOD 

The explicit form of the extended Holstein-extended 
Hubbard (EHEH) model is given by the Hamiltonian 

Heheh = ~* / jCq+i.oA.g + H.c.) +uj 22 a i a i 

i,cr i 

i i 
+ 9opY^H a l + a i) + (4+1 + £li+l)R,<r, (2) 

The notation is the same as in eqn. [TJ In the remainder 
of the paper, we will measure all energies in units of t. 
This model is similar the "extended Holstein- Hubbard" 
model studied in the case of bipolarons^ 3 - 

Working in the path-integral representation and trac- 
ing out the phonon fields, we can express the action of 
the EHEH model as, 

s = I v4( iw - fc)Vv,fc 

Ja.k 



,a',{k} 



g(fa,k2, fc 3 , ki)^, k ^i M ^aM^a' M>(?>) 



where the coupling function g is given by 



g(k 1 ,k 2 ,k Zl k i ) 
2g e 2 p^o 



[ujI + (u>i - w 3 ) 



U + 2Vcos(k 3 
r l + cos(ks - 



-fci) 
fci)] 



(4) 



To analyze this Hamiltonian, we will use our ver- 
sion of the multi-scale functional renormalization group 
(MFRG) method.—^ This method has been used suc- 
cessfully to study a number of problems£&i£ in models 
involving both e-ph and e-e interactions because it cap- 
tures not only the electron correlation effects but also the 
coupling of different energy scales induced by phonon re- 
tardation effects. In its simplest form, the MFRG can be 
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used to follow the flows of the various couplings as the en- 
ergy cut-off is lowered. However, it has been shown that 
retardation can modify the various scalings between dif- 
ferent correlation functions for the Luttinger liquid^ ' 48 ' 49 
Therefore, instead of studying the couplings, we use the 
MFRG to study directly the various susceptibilities which 
take the contributions from different energy scales into 
account XL The dominant ordering can be inferred from 
the flows of the susceptibilities corresponding to different 
order parameters as the cutoff is lowered. 

In the ensuing section, we calculate the susceptibilities 
corresponding to the charge density wave (CDW), the 
spin density wave (SDW), and both the the singlet (SS) 
and triplet (TS) superconducting order parameters. The 
dominant instability among them is determined by the 
most divergent susceptibility as the cut-off, A, is lowered. 
For completeness, we present in the Appendix a short 
summary of the MFRG method and the scaling equations 
for the susceptibilities. For further details, readers are 
referred to several previous articles . 7 ' 17 ! 46 ' 47 In the next 
section, we turn directly to the results of our MFRG 
study. 



III. RESULTS FOR THE SUSCEPTIBILITIES 
FLOWS 

In Fig. [1] we show the susceptibility flows for three 
different sets of e-e interaction parameters, for fixed e-ph 
coupling and phonon frequency. In the top panel of Fig. 
[TJ we choose values of U (=1.00) and V (=0.60) which 
are in the region where, for the extended Hubbard model 
without phonons, the dominant susceptibility is known 
to be a CDWJi Our calculations show that this result is 
unchanged by the inclusion of phonons for the values of 
g ep and ujq shown in the figure. 

In the middle panel, we show the flows for U = 1.00 
and V = 0.20 in the region where, for the extended Hub- 
bard model without phonons, the dominant susceptibility 
is known to be a SDWf^ 7 - Again, we find that this result 
is unchanged by the inclusion of phonons for the values 
of g ep and u>o shown in the figure. 

Finally, in the bottom panel of Fig. Q] we tune down 
the e-e interaction to U = 0.20 and V = 0.20. In the 
extended Hubbard model without phonons, this would be 
in the CDW region. Instead, we find that for the values 
of the phonon parameters in Fig. [TJ the superconducting 
susceptibility clearly becomes dominant. This provides 
an in-principle proof that non-local e-ph interactions are 
able to drive the pairing instability even at half-filling. 

In Fig. [2] we display the "phase diagram" in the U,V 
plane for e-ph coupling, g ep = 0.5 with u>q = 1.0. The 
phase diagrams show that three possible phases-CDW, 
SDW, and SS-are present, with the SS phase being re- 
stricted to small values of U and V, as one would expect. 
For larger values of U and V and for these values of the 
phonon parameters, the model exhibits physics similar to 
that of the conventional extended Hubbard models We 
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FIG. 1. Susceptibility flows for ujq — 1.0, g ep = 0.5 with 
various electron-electron interactions U and V. U = 1.0, V = 
0.6 for the top panel, U — 1.0, V = 0.2 for the middle panel, 
and U — 0.2, V = 0.2 for the bottom panel. 



note that the delicate BOW phase mentioned in the in- 
troduction that occurs between the CDW and the SDW 
in the extended Hubbard model is not shown in our Fig. 
[5] From our previous studies^ 7 - we know that this phase is 
not captured by the MFRG with Nj, — 2 (see Appendix) , 
so we do not expect it to appear in our calculation. This 
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FIG. 2. Phase diagram for g ep — 0.5 with uiq = 1.0. 



choice was made deliberately to allow us to focus on the 
effects of e-ph interactions and phonon retardation, which 
are the phenomena that drive the superconducting fluc- 
tuations. As one would expect intuitively, the most in- 
teresting regime for SC is at weak e-e repulsion. In this 
region, we find the SS susceptibility unambiguously dom- 
inates. This is not what we see in the usual models with 
local e-ph interactions and this arises from the non-local 
nature of the e-ph interactions. 



IV. ANALYTIC INSIGHT INTO THE GAPLESS 
SC PHASE IN THE EHEH MODEL 

In order to gain insight into the mechanism driving 
the dominant superconducting fluctuations in more de- 
tail, we next undertake an essentially analytic study that 
uses a standard renormalization group (RG) approach to 
the bosonized form of the EHEH model. For ease of pre- 
sentation and comparison to the prior literature, we use 
the RG formulation of Yonemitsu and Imada 39 to write 
the EHEH Hamiltonian entirely in terms of bosonic fields. 
Here we will keep only the lowest energy electrons-those 
around the Fermi points at kf and — fc/-and the single 
phonon excitations that couple to the electrons (so the 
2k f and momentum phonons). Because the difference 
between the EHEH model and previously studied models 
lies entirely in the electron-phonon interactions, we will 
focus first on their effects. Explicitly, this means that 
we ignore the back scattering and the Umklapp scatter- 
ing from the direct e-e interactions and incorporate the 
remaining forward scattering into the Luttinger parame- 
ters. The purely electronic part is given by the standard 
free boson Hamiltonian, which we do not need here. With 
these simplifications, the e-ph interaction can be written 



H{ v = / dx— {~f 1 exp[V2i<S> p (x)]cos[V2$ a (x)} 
J na 

4>2k f {x) + H.c.} 
H 2 P = 12 I dx[p + {x) + p^(x)](f> {x), 



Here 71 = (1 + exp(i2k f ))g ep , 72 = 2g ep ; cf> , and <p 2 k f 
are the k = and the k = 2k f components of the phonon 
fields respectively; $ p , and are the bosonized fields for 
the charge part and the spin parts, respectively; p + , and 
P- are the charge densities at the k — kf, and k = —kf 
respectively; and a is the long distance cutoff. Note that 
in the HEH model, 71 = 72 = g e p- At this point, it is 
clear that at half- filling where kf = tt/2 the e-ph inter- 
action contributes only to the linear forward scattering, 
H2. The contribution from the non-linear term, Hi, van- 
ishes exactly at half-filling. This is essentially the same 
result as shown in Eqn. (U). Intuitively, this implies that 
the zero wave vector ordering benefits from the forward 
scattering but the 2k f ordering, like that of the charge 
ordering, does not. Therefore, a charge gapless SC order- 
ing is preferred to the charge gapped CDW ordering. The 
vanishing of Hi at half-filling also means that the charge 
part and the spin parts are decoupled in the absence of 
direct e-e interactions. As a result of these peculiar fea- 
tures, this model can be diagonalized exactly, and it can 
be shown that the charge gap is zero4£ 

When we include direct e-e interactions, the bosonized 
Hamiltonian is no longer solvable. But the special fea- 
ture of this model-namely, the spin-charge separation- is 
preserved. Thus the spin part of the Hamiltonian will not 
be affected by the e-ph interaction which cou ples solely 
to the charge. In addition, as shown in Ref. [39(, since 
the operator (proportional to V) for the spin-charge cou- 
pling in the extended Hubbard model has scaling dimen- 
sion 4, the charge part can be understood independently 
from the spin part for small coupling. Hence to illustrate 
the nature of gapless phase under the influence of both 
e-e and e-ph interactions, we employ the RG equations 
in the form used by Yonemitsu and Imada 39 and focus 
solely on the flows of the charge couplings. As shown in 
Ref. [HH, these coupled flows are described by variables 
X p = 2(1 - k,- 1 ), where k p is the usual Luttinger Liquid 
charge exponent and Y p = g^/^nvf), where 53 is the usual 
Umklapp coupling in g-ology^ From the general form of 
Eqns. (3.11) and (3.12) in Ref. (39|, we consider the 
special case of 1/2 filling and vanishing backward scat- 
tering from the phonons (see Eqn. [S]). In this case, the 
equations in Ref. [391] simplify considerably and become 



dXpjl) 

dl 
dY p {l) 

dl 



-Y p 2 (l) + 2Y 2 (l)D (l), 
-X p (l)Yp(l). 



(6) 



H ep 



Hl p 



rjep 
n 2 1 



(5) 



Here the phonon propagator at momentum is Dq — 
J^yexp(--|a_), with E(l) = E f exp(-l). E f is the 
Fermi energy and Vf is the Fermi velocity. The ini- 
tial conditions are X p (0) = -(U + 6V)/(nvf), Y p (0) = 
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(U - 2V)/(wVf), and F 2 (0) = 2g% p /(m>fU$). Notice that 
Y2 does not flow with They are renormalization group 
equations for the sine-Gordon model, with the exception 
of an additional drift term from the forward scattering 
of the e-ph interaction scaled by the phonon propaga- 
tor. The effect from the phonon propagator exists only 
roughly above the phonon energy. If one momentar- 
ily ignores the e-ph interaction term in the first equa- 
tion, these equations can be integrated exactly to give 
(Xp—Y 2 ) = C, which shows that the RG trajectories are 
hyperbolas in the X p — Y p plane. Further, there is a line 
of fixed points at Y p —0 which are attractive (stable) for 
X p > and repulsive (unstable) for X p < 0. Hence the 
RG flows are as shown in the upper panel of Fig. [3] In the 
quadrant bounded by X p = ±Y p and X p > 0, the flows 
go to a fixed point along the line Y p = 0, which corre- 
sponds to the gapless phase in the sine-Gordon model and 
to the regime of dominant SC fluctuations in the EHEH 
model. Outside this region, the equations always flow to 
strong coupling, reflecting instabilities toward the SDW 
(Y p -> 00, X p -> -co) and CDW (Y p -> -00, X p -> -00) 
phases in the EHEH model . 11 ' 50 In the lower panel of Fig. 
131 the results of including the e-ph coupling are shown. 
Since the RG equation is now explicitly "time" (that is, 
scale) dependent from the phonon propagator, the RG 
flows not only depend on the initial values of the X p and 
Y p , but also on the initial scale, I = Iq, Therefore, we 
cannot draw the RG trajectories as in the case without 
e-ph interaction, and we cannot solve the problem ana- 
lytically. Instead, we use numerics and illustrate some 
typical flows near the boundary between the gapless and 
gapped phases as functions of g cp , with ojq set to 1. The 
effect of the e-ph interaction is to "pull" the flows to- 
wards a larger value of X p . If X p is pulled far enough, it 
will eventually flow to a finite fixed point, while Y p will 
flow to zero (see fig. [3]), and the system will be in the 
region of dominant SC fluctuations. Thus initial values 
that for small e-ph coupling flow to one of the density 
wave phases; for larger e-ph coupling flow instead to the 
superconducting phase. These approximated mostly an- 
alytic results thus corroborate our MFRG study which 
show that a charge gapless phase with dominant super- 
conducting fluctuations can be obtained in the EHEH 
even with all repulsive e-e interactions where both U and 
V are positive. 



V. CONCLUSIONS 

We have shown that by incorporating non-local 
electron-phonon interactions we can find a ID 
Hamiltonian-the EHEH model-that exhibits a gap- 
less metallic phase with dominant superconducting 
fluctuations over a wide range of parameters. This 
contrasts with the results for models involving local 
electron-phonon interactions, where a gapless metallic 
phase, if it exists, does not appear to have dominant 
superconducting fluctuations and occurs over a limited 





FIG. 3. The flows of the charge coupling, X p , and Y p without 
(the upper panel) and with electron-phonon couplings (the 
lower panel). See section IV for a detailed description. 



parameter range. Although a true superconducting 
phase cannot exist in ID, our results suggest that the 
inclusion of non-local e-ph interaction may be crucial to 
understanding the behavior of the quasi-two dimensional 
high-Tc materials 5 -^ and organic superconductors 4^ 
Indeed, the recent experimental evidence that bosonic 
modes, very likely related to phonons, play a crucial role 
in some of the most interesting features of the high T c 
cuprate superconductors 51 provides further motivation 
for studying models with both e-e and non-local e-ph in- 
teractions. In this regard, we point out the recent study 
of the two-dimensional case for c-ph interactions beyond 
the local Holstein coupling where it has been shown 
that superconductivity does exist with some anisotropic 
non-local coupling even at half-filling; moreover, for this 
two-dimensional system, there is very likely to be a true 
long range SC order at zero temperature^ 
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VII. APPENDIX 

A. The MFRG method and equations for the flows 
of susceptibilities 
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tivity (TS) are given by: 

X S A = jD(l,2)f( Pl )f( P2 )(c pui c- Put cl p2t clj. (9) 




The functional renormalization group for electron- 
electron coupled systems has been discussed exten- 
sively in the liteature^2r— The MFRG approach 
implementatio n 46 ! 47 at the one-loop level yields the fol- 
lowing RG flow equations for the coupling functions, 
g(k\,k2,k$,ki), with initial conditions as given by Eqn. 
H: 

dg(ki,k 2 ,k 3 ) _ 
dA 

dp-nr[GA(p)GA(k)]g(k 1 ,k 2 ,k)g(p, k, fc 3 ) 

dp-^[G A {p)G A (qi)]g(p, k 2 ,qi)g(ki,q 1 ,k 3 ) 

dp-^[G A {p)G A (q 2 )][-2g(ki,p, q 2 )g(q 2 ,k 2 ,k 3 ) 
+g{p, fa, 92)5(92, k 2 , k 3 )+g(kx,p, q 2 )g(k 2 ,q 2 , k 3 )}, (7) 



where k = ki + k 2 — p, qi = p + k 3 — ki, q 2 = p + 
k 3 — k 2 , J dp = J dpJ2uj 1 / (27T/3), and G A is the self- 
energy corrected propagator at energy cut-off A. Since 
the interaction vertices are frequency dependent, there 
are also self-energy corrections. At the one-loop level, 
the self-energy MFRG equation is: 



= ~f dp JK [GA {P)] iMP,k,k)-g(k,p,k)}.(8) 



FIG. 4. Discretization of the momenta in the Brillouin zone 
and frequencies in the frequencies axis. This figure shows the 
case N k = 2, N w = 15. 



For 8 = SS, f(p) = 1, whereas for 8 = TS, f(p) = sin{p). 

The static 2k f charge density wave (CDW) and spin 
density wave (SDW) susceptibilities can be written as 7 -: 



Xa= 



fl,CT2 

For 8 = SDW: s t = l,sj. = -1, and for 8 = CDW: 
s t = l,s i = 1. 

In the above equations, Pi is the momentum at energy 
fc, / D(l, 2) = Jj fi|>A d£i J(6) 4 2|>A dfc J(6), and J(£ 
is the Jacobian for the coordinate transformation from k 
to £fc. The dominant instability is determined by the 
most divergent susceptibility as the cut-off A is lowered. 

The RG flows for the SS and TS susceptibilities are, 



dA 
dA 



'dp—[G A {p)G A {-p)]{z{( P )Y, (11) 



W^\GA{p')G A {-p')]Zi{p') g y,p), (12) 



where g s (p',p) = g(p', -p', -p) for both 8 =SS and TS. 
The RG flows for the SDW and CDW susceptibilities 



We have solved the coupled integral-differential equa- 
tions, (JTJ) and ©, numerically with two Fermi points 
(Nk = 2) and by dividing the frequency axis into fifteen 
segments (N u — 15). Fig. @] shows the discretization 
scheme for Nk — 2 and N u = 15. 

We next calculate within the MFRG approach for the 
RG flow of susceptibilities in the static (zero frequency) 
and long- wavelength limit. The static susceptibilities for 
singlet superconductivity (SS) and triplet superconduc- 



dA 



' dp dA [GA(p)GA(p+Q)]( ^ (p))2 ' (13) 



^ W -=/dj/^[G A (y/)r; A (//+^,]Z;i(y/lr/(y/. y ,)il4l 



dA 



where Q = (Q = vr,0). For 8 = SDW: g s (p',p) = ~g{p + 
Q,p',p), and for 8 = CDW : g s (p',p) = 2g(p',p + Q,p) - 
g(p + Q,p' ,p). The function Z s (p) is the effective vertex 
in the definition of the susceptibility \ S ^ 
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